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I— I Abstract 

o 

\^ The symmetron is a scalar field associated with the dark sector whose coupling to mat- 

i-^ ter depends on the ambient matter density. The symmetron is decoupled and screened in 

JL regions of high density, thereby satisfying local constraints from tests of gravity, but couples 

^ with gravitational strength in regions of low density, such as the cosmos. In this paper we 

d derive the cosmological expansion history in the presence of a symmetron field, tracking the 

evolution through the inflationary, radiation- and matter-dominated epochs, using a combi- 

^ nation of analytical approximations and numerical integration. For a broad range of initial 

,__! conditions at the onset of inflation, the scalar field reaches its symmetry-breaking vacuum 

by the present epoch, as assumed in the local analysis of spherically-symmetric solutions 



(N 

I — ^ and tests of gravity. For the simplest form of the potential, the energy scale is too small for 

^■-^ the symmetron to act as dark energy, hence we must add a cosmological constant to drive 

^~~* late-time cosmic acceleration. We briefly discuss a class of generalized, non-renormalizable 

. 5^ potentials that can have a greater impact on the late-time cosmology, though cosmic accel- 

rN eration requires a delicate tuning of parameters in this case. 



1 Introduction 

The standard ACDM cosmology explains cosmic acceleration with a cosmological constant, a 
logically sound but seemingly highly tuned explanation. Therefore there has been much interest 
in exploring the possibility that the cosmic acceleration could be caused by a previously unobserved 
dynamical component of the universe [1] — [5]. At the same time, nearly massless gravitationally 
coupled scalars are generically predicted to exist by many theories of high energy physics. No 
experimental sign of such a fundamental scalar has yet been seen, in spite of tests designed to 
detect solar system effects and fifth forces that would generally be expected if such scalars were 
present [6, 7]. 

If light scalars exist, they must utilize a screening mechanism [8] to hide themselves from 
local experiments. Screening mechanisms rely on non-linearities whose behavior depends on the 
ambient matter density. In regions of high density, such as the local environment, the scalars 
develop non-linearities that effectively decouple them from matter. In regions of low density, such 
as the cosmos, the scalars couple to matter with gravitational strength and mediate a long-range 
force, thereby affecting the nature of gravity and the growth of structure on large scales [9]. 

Aside from cosmology, these theories find independent motivation in the vast experimental 
effort aimed at testing the fundamental nature of gravity at long wavelengths [7]. Viable screening 
theories make novel predictions for local gravitational experiments. The subtle nature of these 
signals have forced experimentalists to rethink the implications of their data [10, 11] and have 
inspired the design of novel experimental [12] — [16] and observational tests [17] — [21]. The theories 
of interest offer a rich spectrum of testable predictions, from laboratory to extra-galactic scales. 

Only a handful of successful and robust screening mechanisms have been proposed to date [8] . 
The first is the Vainshtein mechanism [34, 35, 36], which works when the scalars have derivative 
self-couplings which become important in the vicinity of massive sources such as the Earth. The 
strong coupling boosts the kinetic terms, so that after canonical normalization the coupling of 
fluctuations to matter is weakened. This mechanism is central to the phenomenological viability 
of brane- world modifications of gravity [37] — [52], massive gravity [36], [53] — [61] (see [62] for a 
review), degravitation models [63] — [66] and galileon theories [67] — [80]. See [81] — [86] for some 
phenomenological implications. 

A second screening mechanism is the chameleon mechanism [22] — [31], which works when the 
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scalars are non-minimally coupled to matter in such a way that their effective mass depends on 
the local matter density. Deep in space, where the local mass density is low, the scalars are light 
and display their effects, but near the Earth, where experiments are performed, and where the 
local mass density is high, they acquire a mass, making their effects short range and unobservable. 
The effective coarse-grained description of chameleon theories, including careful considerations of 
averaging, has been derived in [28]. Chameleonic vector fields, such as gauged B — L, have been 
proposed in [32]. See [33] for an attempt at realizing the chameleon in string compactifications. 

Recently, two of us have proposed a third screening mechanism, called the symmetron [87], 
based in part on earlier work [88, 89]. For this mechanism to operate, the vacuum expectation 
value (VEV) of the scalar must depend on the local mass density. The VEV becomes large in 
regions of low mass density, and small in regions of high mass density. In addition, the coupling 
of the scalar to matter is proportional to the VEV, so that the scalar couples with gravitational 
strength in regions of low density, but is decoupled and screened in regions of high density. 

This is achieved through the interplay of a symmetry-breaking potential [87], 

n0) = -^/^v+^A0^ (1) 

and a Z2-invariant universal conformal coupling to the trace of the matter stress tensor, ~ 
(jP'T^ JM"^, so a local matter density contributes to the effective mass of the scalar. In vacuum, 
the scalar acquires a VEV |0| = 0o = /w/vA, which spontaneously breaks the reflection symmetry 
of the lagrangian — )■ —0. In the presence of sufficiently high ambient density, on the other 
hand, the potential does not break the symmetry and the scalar is trapped near = 0. This 
is shown in Figure 1. In addition, 50 fluctuations couple to matter as (0vev/^^)<50 P- Hence, 
symmetron perturbations are weakly coupled in high density backgrounds and relatively more 
strongly coupled in low density backgrounds. 

The symmetron naturally takes the form of an effective field theory. The potential comprises 
the most general renormalizable terms invariant under the Z2 symmetry — )■ —0. The coupling 
to matter is the leading such coupling compatible with the symmetry. It is non-renormalizable, 
suppressed by the mass scale M, thus the symmetron is an effective theory with cutoff M . Re- 
markably, tests of gravity constrain this cutoff to be around the GUT scale [87], so any GUT 
theory with a low energy scalar might be expected to yield a symmetron-type lagrangian at low 
energies. From this point of view, symmetron theories are more natural-looking than chameleon 

2 



models. As with all non-supersymmetric scalars, however, the coupling to matter generates large 
quantum corrections to the mass which must be fine-tuned away. 

Symmetron theories predict a host of observational predictions, some of which are distin- 
guishable from other screening mechanisms. Although the local environment is dense, and there- 
fore one in which the symmetron-mediated force is weak, the symmetron nevertheless leads to 
small deviations from General Relativity in the solar system. As reviewed in Section 2.3, the 
predicted signals for Lunar Laser Ranging and Mercury's perihelion precession are just below 
current bounds and within reach of next-generation experiments [87]. On the other hand, the 
signal from binary pulsars is much weaker and should not be observed, since both neutron stars 
are screened. This makes the symmetron distinguishable from over-the-counter scalar-tensor the- 
ories, such as Brans-Dicke (BD) theories, for which solar system and binary pulsar signals are 
comparable. The symmetron solar system predictions are also distinguishable from chameleon 
and Vainshtein screening [87]. 

SjTiimetrons, like chameleons, predict effective macroscopic violations of the equivalence 
principle between large (screened) and small (unscreened) bodies. Dwarf galaxies in low-density 
environments, in particular, offer a fertile playground for testing these ideas [18, 21]. The stars are 
oblivious to the symmetron whereas the hydrogen gas experiences an additional force, resulting 
in an enhanced rotational velocity for the gas [18]. Meanwhile, the infall motion of a dwarf galaxy 
can lead to a segregation of the stellar disk from the dark matter and the hydrogen gas [21]. 

We begin with an extensive overview of symmetron physics in Section 2, reviewing the 
spherically-symmetric solution (Section 2.2) and constraints from solar system tests of gravity 
(Section 2.3). Section 2.1, in particular, offers a novel perspective of the symmetron in terms 
of Jordan-frame variables. In this frame, the action describes a Brans-Dicke type scalar-tensor 
theory, with a field-dependent Brans-Dicke parameter. The screening mechanism is understood 
as the Brans-Dicke parameter becoming large in regions of high-density, corresponding to the 
decoupling of the scalar. In this respect, the symmetron is qualitatively similar to the Vainshtein 
mechanism, though its non-linearities stem from potentials rather than derivative self-couplings. 

The remainder of the paper focuses on the cosmological evolution in the presence of a sym- 
metron field. Specifically, we are interested in checking whether a symmetron leads to an ac- 
ceptable cosmology and, more interestingly, whether it can act as dark energy driving cosmic 
acceleration. We derive the evolution of the scalar field from the infiationary epoch until the 
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(a) Effective potential in regions of high am- 
bient density. The VEV is 0. 



(b) Effective potential in regions of low am- 
bient density. 



Figure 1: Schematic plots of the syinmetron effective potential, illustrating the symmetry breaking 
phase transition. 



present. The dependence of the symmetron effective potential on the matter density has implica- 
tions for the cosmological evolution of the scalar field. Since the matter density redshifts in time, 
the effective potential is time-dependent and results in a phase transition when the matter density 
falls below a critical value. We will choose parameters such that the phase transition occurs in 
the recent past (-Ztran ^ !)• 

A key question is whether the evolution allows the scalar field to reach the symmetry- 
breaking vacuum by the present epoch, as assumed in the analysis of solar system tests and other 
phenomenological studies. The short answer is that the coupling to matter efficiently drives the 
field towards the symmetry-restoring point, so that for a broad range of initial conditions it reaches 
(f) c^ well before the phase transition. How the field makes it there, however, is an interesting 
story, as summarized below. 

In Section 3 we describe the symmetron evolution during the radiation- and matter-dominated 
eras of standard big bang cosmology. The effective mass squared of the symmetron due to its cou- 
pling to matter fields is m^g ~ T^/M'^ ~ (1 — 3w)p/M'^, where w is the equation of state. During 
the radiation-dominated era, T^ is dominated by the (subdominant) non-relativistic component, 
hence rn^Q ~ Q^H'^Mp^/M'^. At early times, therefore, m^Q <^ H^, and the symmetron remains 
essentially frozen at its initial value, denoted by 0rad-i- As the universe cools and flrn increases, 
eventually m^g ~ H^, and the symmetron starts rolling and undergoes damped oscillations around 
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= 0. Since M <^ Mpi, this rolling phase is triggered deep in the radiation-dominated epoch, prior 
to Big Bang Nucleosynthesis (BBN). From this moment onwards, the field amplitude decreases by 
a factor of (M/Mpi)^/^ until matter-radiation equality. The damped oscillations continue during 
the matter-dominated era, and the field amplitude decreases by an additional factor of 10"'^ until 
the present epoch. Thus, for a broad range of initial field values, the symmetron ends up close to 
the origin by the onset of the phase transition. 

In Section 4, we follow the evolution of the symmetron through the phase transition. For 
the quartic potential (1), the symmetron nicely tracks the effective minimum and settles to the 
symmetry-breaking vacuum by the present time. However, the energy scale of the potential is too 
small to drive cosmic acceleration, and the symmetron backreaction is negligible. Hence we must 
add a suitable cosmological constant, thereby making the expansion history indistinguishable from 
ACDM cosmology. 

We also consider generalized, non-renormalizable forms of the potential, with scales chosen 
to have a greater impact on late-time cosmology: i) the potential energy difference between 
maximum and minima is of order HqMp^, as desired to act as dark energy; and ii) the mass of 
small fluctuations around the minima is of order Hq, such that the symmetron can impact the 
growth of structure on the largest observable scales. The resulting potential is sharply peaked 
at = 0, and displays very shallow minima. Because of this asymmetric form, shortly after the 
onset of the phase transition the symmetron overshoots the minimum and reaches a value of order 
Mpi. As a result, the fleld does not converge to the symmetry-breaking vacuum by the present 
time, which is problematic for tests of gravity. Meanwhile, the symmetron energy density quickly 
becomes irrelevant after the transition. One way to circumvent these problems is by requiring the 
phase transition to occur in the very recent past (ztran ^ 1); but this requires a delicate choice of 
parameters. 

In Section 5 we go back to the early universe and describe the evolution of the symmetron dur- 
ing a period of inflation and subsequent reheating. Since the effective mass mlg ~ 12H^^^Mpi/M'^ 
is ^ if^f during inflation, the scalar undergoes rapid oscillations with exponentially decaying 
amplitude, thanks to Hubble damping. The end result is that a fleld initially displaced from 
= ends up exponentially close to this point. At the end of inflation, however, as the inflaton 
oscillates around its minimum, its equation of state also oscillates between w = —1 and w = 1. 
In turn, ml^ oscillates between negative and positive values, which can result in ampliflcation. 



With a simple parainetrization of w{t), we find that the symmetron evolution can be de- 
scribed at early times by a Mathieu equation, with solution given by the Mathieu cosine function. 
As usual with Mathieu functions, the stability behavior depends on parameter choices, in this case 
the inflation mass (setting the rate of oscillations in wi^g), its decay rate (setting the duration of 
the unstable phase) and the effective symmetron mass during inflation. We place a constraint on 
these parameters by demanding that the symmetron displacement be within effective field theory 
bounds. We summarize our results and discuss future directions in Section 6. 

2 Review of the Symmetron 

The symmetron is a special case of a general scalar-tensor theory, 

/r/\/f2 1 1 /■ , 

d'x V^ -^R - -g'-'d.^d,^ - V{^) + J d'x v^An iij, ~g,,) , (2) 

where the metric has mostly-plus signature, and £m is the lagrangian for the matter fields, denoted 
collectively by ip. The matter fields couple minimally to the Jordan frame metric, 

^^^ = A^(0)^^^, (3) 

related to the Einstein frame metric g^^, by a positive function A{(j)). In particular, (j) couples 
universally to all matter fields, hence the weak equivalence principle holds. See [90] for a discussion 
of quantum corrections to universal coupling. 

The scalar field equation of motion is given by 

n0 = y<^-A3(0)A^(0)T, (4) 

where T = g^^T^^ is the trace of the the matter stress-energy tensor, T^^, = —{2/y/^)5C^/5g^^, 
covariantly conserved with respect to g^j^y-. V^T^^^ = 0. Meanwhile, the Einstein equations are 

M^,G,y = T^, + A'i<l>)f,y. (5) 

Note that because (p couples conformally to matter, its stress-energy tensor, 

,,u = c)f,(p(Ju(p - -> 

is of course not covariantly conserved: V^T/f^, ^ 0. 
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Tt = d,<pd,<P - ^9Ad4>Y - g,.V{ct>) , (6) 



The form of the functions A{(j)) and V{(f)) is crucial to the operation of the symmetron 
mechanism. These functions are assumed symmetric under — ;■ —0 and are such that the effective 
symmetry breaking potential (41) has a zero VEV for large p and a large VEV for small p. In 
addition, the function A{(j)) should be such that the coupling of scalar fluctuations to matter is 
proportional to the VEV. 

The simplest symmetron theory, considered in [87], is that of a quartic potential (1) and 
quadratic coupling: 

The potential V{(j)) comprises the most general renormalizable form invariant under the Z2 sym- 
metry (j) — )■ —(f). The coupling to matter ~ cf? /M"^ is the leading such coupling compatible with 
the symmetry. However, as we will see shortly, once we impose constraints from tests of gravity, 
the energy scale in this potential is too small to act as dark energy. In Section 4.2, we will consider 
generalizations that have a greater impact on cosmic acceleration. 

The model involves two mass scales, p and M, and one positive dimensionless coupling A. 
The mass term is tachyonic, so that the Z2 symmetry (j) — ;■ —0 is spontaneously broken. With 
this canonical choice, we will see shortly that the relevant field range is ^ M, such that any 
0(0^ /M'^) terms in v4(0) can be consistently neglected. For the case of non-relativistic matter, 
which is relevant for most applications including local tests of gravity, T ~ —p. Expressing (4) 
in terms of p = 74^(0)p, which is 0-independent, the effective potential is, up to an irrelevant 
constant, 

M0) = ^ (^ - /i') 0' + ^A0^ (8) 

Whether the quadratic term is negative or not, and hence whether the Z2 symmetry is sponta- 
neously broken or not, depends on the local matter density. 

The screening mechanism works roughly as follows: in vacuum or in large voids, where 
p ~ 0, the potential breaks reflection symmetry spontaneously, and the scalar acquires a VEV 
101 = 00 = p/\/~X] in regions of high density, such that p > M'^p'^, the effective potential no longer 
breaks the symmetry, and the VEV goes to zero. Meanwhile, to lowest order the symmetron- 
matter coupling is ~ pcjp' /M"^. Fluctuations 50 around the local background value 0vev, which 
would be detected by local experiments, couple to density as 

^^^^^P- (9) 
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In particular, the coupling is proportional to the local VEV. In high-density environments where 
the symmetry is restored, the VEV should be near zero and fluctuations of do not couple to 
matter. In less dense environments, where p < M'^p? and the symmetry is broken, the coupling 
turns on. 

We will be interested in the case where the field becomes tachyonic around current cosmic 
density, 

HlMl, ~ p^M^ , (10) 

so that the phase transition occurs around the onset of cosmic acceleration. This fixes /i in terms 
of M, and hence the mass itiq of small fluctuations around the symmetry-breaking vacuum: 

mo = 72/. ~ ^H, . (11) 

As reviewed in Section 2.3, constraints from local tests of gravity require M < 10~^Mpi. Hence the 
range m^^ of the symmetron-mediated force in regions of low mass density is <0.1 Mpc — too 
heavy to act as a slowly-rolling quintessence field, but light enough to impact structure formation 
and have interesting astrophysical implications. 

The symmetron-mediated force F^ relative to the Newtonian force F^ between two test 
masses in vacuum is set by the symmetry-breaking value 0o- 
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.. Hl^V. (12) 



If the scalar-mediated force is to be comparable to gravity in vacuum, then we must impose 
0o/M2 ~ 1/Mpi, that is, 

where g ~ 0{1). To be precise, it follows from (12) that g measures the strength of scalar force 
in vacuum relative to gravity: F^ = 2g^F^. 

Combined with (11) and the requirement M < 10~^Mpi, (13) fixes the dimensionless quartic 
coupling to be exponentially small: 

Note that since M ^ Mpi from tests of gravity considerations, (13) implies that <^ M, hence 
the entire field range of interest lies within the regime of the effective field theory. 



Unfortunately, these requirements imply that the potential energy is too small to have sig- 
nificant backreaction and drive cosmic acceleration. Indeed, using (10) and (13), the potential 
height difference between = and = 0o is 



^V 



4A 



HlM^ < HlMl^ . 



In Section 4.2 we will consider generalized forms of the potential that have AV 
therefore better candidates for dark energy apphcations. 



(15) 



HqMp^ and are 



2.1 Jordan-Frame Description 

Because the scalar is assumed to couple universally to matter, we can also understand the sym- 
metron mechanism in Jordan frame. We will see that the coefficient of the symmetron kinetic 
term depends on the local value of the scalar field. This coefficient becomes large in regions of 
high density, which effectively decouples the scalar field from matter. From this point of view, 
the symmetron mechanism is qualitatively similar to Vainshtein screening, though it relies on a 
potential rather than derivative self couplings. 

In terms of the Jordan-frame metric ^^jy = A^ {(f)) g f^,, , the action (2) takes the form 
"M2 



S 



I 



d X W—g 



^^R - ^^^^^r^.'^d,'^ - ^ V(M/)1 + /" d 



X V-9^m{^,9f,u) , (16) 



where "^ = A ^(0). Thus the symmetron behaves as a Brans-Dicke scalar field, with field- 
dependent Brans-Dicke parameter 



^(0) 



1 



2M2j(dlnA/d0)^ 



M2 



Ml 



pi<, 



-3 



(17) 



where in the last step we have substituted the quadratic form A{(j)) ~ 1 + 0^/2M^ and used the 
fact that A{(j)) c^ 1 for all field values of interest. 

The nature of the screening mechanism in this frame is crystal clear. In regions of high 
density, where ~ 0, the Brans-Dicke parameter is large. 



tu 



M' 



2Mi 



pi<; 



>l 



(high density regions) 
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indicating that the scalar field decouples and General Relativity is recovered. In regions of low 
density, where |0| ^ 0o = gM'^/Mpi, the Brans-Dicke parameter is of order unity, 

2 + 3u; ^ — - (low density regions) , (19) 

hence the scalar field generates order unity corrections to General Relativity. 



2.2 Spherically Symmetric Solution 

To study the implications for tests of gravity, we are interested in the symmetron profile around 
astrophysical sources, such as the Sun. For this purpose, we can safely work in the Newtonian 
limit, ignoring non-linear effects in gravity and the back-reaction of the scalar field on the metric. 
Moreover, the source is approximated as static, spherically symmetric and pressureless (T ~ — p). 
It has radius R and homogeneous mass density p, such that p ^ p?M'^. For simplicity, we further 
assume the object lies in vacuum. 

Written in terms of the density p = A^p, which is conserved in Einstein frame, the scalar 
field equation (4) in spherical coordinates reduces to 

^'^+^^'^ = K^ + A0P- (20) 

Analogously to what was done in [22], this radial field equation can be thought of as a fictional 
particle rolling in an inverted effective potential —Vcs{4>), subject to a "friction" term {2/r)d(f)/dr. 
The boundary conditions are that the solution be smooth at the origin, and approach its symmetry- 
breaking value at infinity: 

^0(0) = 0; 0(r^oo) = 0o. (21) 

First consider the solution inside the object. Since the objects of interest are much denser 
than the current critical density, p ^ p^M"^ ~ H^Mp^, the effective potential can be approximated 
as Kff(0) — pcf? /2M'^. The interior solution satisfying the first of (21) is 

(/)in(r) ~ C - sinh i^r\ (for r < R) . (22) 

This involves one undetermined constant C to be fixed shortly. 

Next consider the exterior solution. It turns out that the potential is negligible, and the 
symmetron evolves as a free field, until it reaches the vicinity of the minimum. In the mechanical 
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analogy, this corresponds to the fictional particle having a large velocity as it exits the object and 
begins its climb towards the maximum of the inverted potential. Since the potential is irrelevant 



except near = 0o, we can make the quadratic approximation \4flf(0) 
r > R. The solution satisfying the second of (21) is then given by 



mn - 



R 

r 



D -e-^'oir-R) ^ 0^ 



(for r > R) 



72 for all 



(23) 



with one undetermined constant D. 

The coefficients C and D are fixed by matching the field and its radial derivatives at the 
interface r = R. In doing the matching, we can safely assume that moR <^ 1, since the symmetron 
Compton wavelength uiq^ is cosmologically large (<0.1 Mpc). It is convenient to express the 
result as 



C 



D 



'AR , 
^OA/ — sech 

1- 



R 
AR 



AR 



tanh 



R 
AR 



(24) 



where we have introduced the thin-shell factor, 

AR M2 M2 



R pR^ 6M2j$ 6^Mpi<l> ■ 



(25) 



Here $ = pR'^/GMp^ denoting the gravitational potential of the source, and where in the last step 
we have used (13). As the notation suggests, AR/ R will soon be interpreted as a thin shell factor 
for the solutions, in analogy with the chameleon mechanism. In fact, (25) precisely matches the 
chameleon thin shell expression, with (J)q being interpreted as the ambient chameleon value. Hence 
one can immediately anticipate that symmetrons and chameleons will have similar phenomenology, 
in particular for astrophysical tests [18, 21]. 

Consider a test particle a distance R <^ r <^ rriQ^ away from the object, such that ^ 0o 
and d\nA/d(j) ~ g/Mpi. Substituting the expression for D into (23) and using (25), the scalar 
force on this particle relative to gravity is 



Fa 



g d0/dr 



Ml 



Qg' 



PI 



, AR 



AR 



tanh 



A 
A^ 



(26) 
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Aside from the coefficient (7, the magnitude of the scalar force is thus determined by a single 
parameter, AR/R. 

Screened objects have relatively large gravitational potential such that AR/R ^ 1. In this 
case the tanh term can be neglected and (26) reduces to 



^ 
Fn 



^6<72^«1. (27) 

screened 



The scalar force between a screened object and a test particle is therefore suppressed by a thin-shell 
factor AR/R ^ 1. This can be understood as an analogue of the thin-shell effect of chameleon 
models [22, 23]. It is clear from the interior profile, 

■ sinh 



^ ^ cosh 



that the field is exponentially suppressed compared to (po everywhere inside the object, except 
within a thin-shell AR beneath the surface. The symmetron is thus weakly coupled to the core 
of the object, hence its exterior profile is dominated by the thin shell contribution. 

Unscreened objects, on the other hand, have relatively small gravitational potential such that 
AR/R in (25) is formally ^ 1. In this regime, we can Taylor expand the tanh and (26) simplifies 
to 

^ ^ 2/ . (29) 

There is no thin shell in this case, and the result is consistent with the scalar force between two 
test masses given by (12) and (13). Indeed, we see from (28) that ~ 0o everywhere inside the 
object, hence the symmetron couples with gravitational strength to the entire source. 

2.3 Constraints from Tests of Gravity 

For the theory to be phenomenologically viable, it must of course satisfy all constraints from tests 
of gravity. Here we review and expand on the original analysis of [87]. Because the symmetron has 
a long Compton wavelength in all situations of interest, tests of the inverse-square-law are trivially 
satisfied. And because its coupling to matter is universal, the weak equivalence principle is also 
satisfied. In fact, as shown in Section 2.1, the symmetron can be understood as a Brans-Dicke 
scalar field with an effectively density-dependent Brans-Dicke parameter. Hence the relevant tests 
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are the same that apply to standard Brans-Dicke theories, namely post-Newtonian tests in the 
solar system and binary pulsar observations. 

What matters for solar system tests is the local field value, since this determines the coupling 
of the symmetron to matter. At a generic point in the solar system, this is determined by the 
symmetron profile interior to the galaxy. Clearly a necessary condition to satisfy local tests is that 
the Milky Way galaxy be screened, ARq/Rq <^ 1, in which case the local value is given by (28) 



Rg Rus cosh^ ^^^^ 



Rg 

where Rg ~ 100 kpc is the Milky Way radius, and R^s ~ 10 kpc is our distance from the galactic 
center. Recall that for the quartic potential considered here, 0o = gM'^/Mp\. We will see that 
post-Newtonian constraints in the solar system translate to a bound on ARq/Rq, and hence on 
M. 

General scalar-tensor theories of the form (16) have two non-vanishing post-Newtonian pa- 
rameters, f3 and 7, defined in terms of the Jordan-frame metric as 

^00 = - (1 + 2$j + 2/3$^) ; 

~gij = (l-2$j7)5,,. (31) 

In terms of uj{(f)), these parameters are given by [7] 

1+W(0) 
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2 + w(0) 



/3 = 1 + 1 -TTY^- (32) 

(3 + 2a;)' (4 + 2a;) d(A-2) ^ ^ 



From (18), the local value of these parameters satisfy 

2Mpi</)g\ 2 /'/'g^ 






The tightest constraint on 7 comes from time-delay measurements with the Cassini spacecraft [91]: 
I7 — 1| < 10~^. The most stringent bound on /3 comes from the Nordvedt effect, which describes 
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the difference in free-fall acceleration of the Moon and the Earth towards the Sun due to scalar- 
induced differences in their gravitational binding energy. Lunar Laser Ranging observations give 
the constraint \f3 -1\< 10''^ [7]. 

Although the experimental bound on /3 is weaker, it actually yields the tightest constraint on 
our model parameters, since the predicted deviation in (33) is ~ Mpj/M^ larger than 7. Assuming 
g ~ 0{1) and using the fact that $g — 10^® for the Milky Way galaxy, we find that |/3 — 1| < 10~^ 
is satisfied for 

^<6xl0-. (34) 

Kg 

In terms of the coupling mass scale M, this translates to 

M < 10"^ Mpi . (35) 

In turn, this corresponds to I7 — 1|g ~ 10^^°, which lies well below the current limit. Intriguingly, 
the upper bound on M is near the GUT scale. Note that if M nearly saturates this bound, then 
objects with $ > 10^®, such as the Sun ($0 ~ 10^^) will be screened but those with $ < 10^^, 
such as the Earth ($9 ~ 10~^), will not. 

If the symmetron is a low energy scalar associated with GUT-scale physics, such that M 
nearly saturates (35), then near- future experiments probing post-Newtonian corrections can poten- 
tially detect symmetron effects. The APOLLO observatory, in particular, should improve Lunar 
Laser Ranging constraints by an order of magnitude [92]. As we have seen, the predicted signal 
for time-delay and light-deflection are much more suppressed. 

Meanwhile, the constraints from binary pulsars are trivially satisfied, since both the neutron 
star and its companion are screened. The force between these bodies is therefore suppressed by 
two thin-shell factors: 

F^ fAR\ fAR\ .fARcY in . ^ 

xK^ -10 Kt^ <10"'°, (36) 



-f^N \ ^ / pulsar V -''' / companion V """G 

where we have estimated $puisar ~ 0.1 ~ IO^'I'g and "^companion ~ 10~^ ~ $0? and the last step 
follows from (34). This is well below the current pulsar constraints on Brans-Dicke scalar-tensor 
theories, which translate to F^/F^ < 10~^ [7]. 
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3 Cosmological Evolution 

Central to the phenomenological viability of symmetron theories is whether the resulting expansion 
and growth histories are consistent with observations. In the remainder of the paper, we will study 
the cosmological evolution of the symmetron field, from the inflationary epoch until today. In this 
Section, we begin by setting up the equations and solving them for the standard radiation and 
matter dominated eras. Given the scales involved, the scalar potential is completely negligible at 
early times, and the symmetron evolution is determined by the interplay of Hubble friction and 
its direct coupling to matter. 

A key question is whether the evolution allows the scalar field to reach the symmetry-breaking 
vacuum by the present epoch, as assumed in the analysis of local tests discussed in Section 2.3. 
We will see that the coupling to matter is efficient in bringing the field close to = by the onset 
of the phase transition (z ~ 1), for a broad range of initial conditions. 

3.1 Preliminaries 

Consider the homogeneous evolution of the scalar field, = 0(t), in a Friedmann- Robertson- 
Walker (FRW) background: ds| = — dt^ + a^(t)dx^. We assume that infiation takes place at 
early times, so that the universe is nearly spatially-flat. The above line element is in Einstein 
frame, hence the label "E". Throughout the analysis, dots will represent derivatives with respect 
to (Einstein-frame) cosmic time t. 

The matter content is modeled as a set of non-interacting perfect fluids indexed by i, each 
with constant equation of state Wi. The trace of the stress tensor therefore becomes T = '^i{—^ + 
3wi)pi. Substituting this into (4), the symmetron equation of motion reduces to 

4> + 3H<p + V4, + A3(0)A,<^ ^(1 - 3wi)pi = . (37) 

i 

By assumption, each matter component is separately conserved: 

p^ ^ a-'^'^""'^ = ^-3(i+-,)a-3(i+«,.) ^ (3g) 

where aj = A{(f))a is the Jordan-frame scale factor. To make the 0-dependence explicit in (37), it 
is convenient to deflne a rescaled energy density, 

^. ^ ^3(i+«,.) ~^^ (39) 
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such that Pi ~ a^'^'^^^'^^'^ satisfies the usual conservation law in Einstein frame. Note that pi does 
not represent the physical density in Einstein frame — it is a mathematical construct, which can 
be treated as independent of cf) when integrating (37). 

In terms of pi, the symmetron equation of motion becomes 

+ 3/70 + V;^ + ^ _^(1 _ 3^.)p^ = . (40) 

i 

For regions of nearly constant pj, we can think of this as a scalar field moving in an effective 
potential 

i 

As the universe expands and the densities pi redshift, this effective potential correspondingly 
evolves in time. Note that the difference between pi and pi is negligible for most of the history of 
the universe, since A{(f)) ^ 1 until the phase transition. 

Meanwhile, the Einstein-frame Friedmann equation, obtained from the (0, 0) component 
of (5), is given by 

?>Ml,H^ = ^02 + Vi^) + A\<P) J2 P^ ■ (42) 

i 

In terms of the 0- independent density pi, this becomes 

3Mi,H' = ^02 + V{<P) + J2 ^'-''"'i^)P^ = I'P' + Kff(0) . (43) 

i 

In summary, the independent cosmological equations are the scalar field equation (40), the Fried- 
mann equation (43) and the conservation laws pi ~ a"^*-^^"''-'. In the remainder of this Section, we 
will study the evolution of the scalar field during the standard radiation- and matter-dominated 
eras. The most recent history, including the phase transition and the onset of cosmic acceleration 
will be discussed in Section 4. 

3.2 Radiation-Dominated Era 

We imagine that after infiation the symmetron ends up at rest at some field value 0rad~i- Later on 
in Section 5 we will estimate 0rad-i by studying the inflationary and reheating dynamics in detail, 
but for the time being let us consider it as a free parameter. To simplify the analysis, we assume 
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that 0rad-i ^ M, SO that higher-order corrections in (7) for the couphng function are neghgible: 

Ai<P) ^ 1 + ^ • (44) 

Moreover, the self-interaction potential V^(0) can be neglected until the late-time phase transition. 
Indeed, given (11) and (14) the quartic potential satisfies \^(0) ~ V{Q) ~ ^i^/^/X < H^Mp^. The 
symmetron kinetic energy is also negligible, as argued shortly. Hence, for most of the evolution 
the symmetron dynamics are governed by a tug-of-war between the coupling to matter driving 
the field and Hubble friction slowing it down. 

Since ^4(0) ^ 1, the Friedmann equation (43) in the radiation-dominated epoch takes the 
standard form 

3H^M^^ = p,{a) , (45) 

where p^ ~ a"'^. Meanwhile, since the radiation stress energy has negligible trace, it drops out of 
the scalar equation of motion, and (40) reduces to 

4> + 3H^ + p^{a)^ = 0, (46) 

where pm ~ a~^. Note that we have neglected the potential term V^, by the argument below (44). 
The physics of this equation is clear — it describes a simple harmonic oscillator subject to friction 
with a weakening spring constant (weakening since pm redshifts in time). At very early times, 
when H^ ^ p^/M^, the oscillator is overdamped, and the field remains stuck at ^ 0rad-i- 
At late times, when H^ ^ Pt^/M"^, the oscillator is underdamped, and the symmetron oscillates 
around its minimum. In other words, we expect the field to remain frozen at its initial value until 
a = aunfreeze, where ttunfreeze Satisfies 

-. PrayOj) <J-'"pi C^unfreeze /'/17^ 

Here Oeq is the scale factor at matter-radiation equality. That is, 

'^unfreeze '^eq „ f,r2 ' \ "J 

Since M^ ^ ^pi) the scalar field begins to roll well before matter-radiation equality. For instance, 
with M = 10""^ Mpi, this occurs at z ~ 10-*^^, i.e., well before BBN. 
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These considerations are borne out by the exphcit solution. Defining 

2y3Mp- / - ^ ^/2 



-^^(^1 ■ («) 



the scalar equation of motion (46) simplifies to 



s -— r + 3s-— + s = 0. (50) 

ds^ as 

The solution with — )• 0rad-i for s ^ 1 is given by a Bessel function: 

0=^^^^Ji(s). (51) 

s 

As advocated, the field remains essentially constant until s ~ 1, which indeed corresponds to 
a ~ ctunfreeze- Subsequently, the field undergoes damped oscillations. At late times s ^ 1, in 
particular, (51) implies 

coss , f M Y^^ ^acq\3/4 f2y/3Mpi f a^ ^^^^ 



^-^-\m;J [f) "M^M^U ■ <^^' 



up to an irrelevant phase. Thus, by matter-radiation equality the amplitude of the scalar field has 
dropped by a factor of 

With M = 10~^Mpi, for instance, this corresponds to a suppression factor of 10~^. 

We are now in a position to check that the symmetron kinetic energy is indeed a negligible 
contribution in the Friedmann equation, as assumed throughout. From the properties of Bessel 
functions, and since H ~ a~^ in the radiation-dominated epoch, we have 



H 2 



d0 



ds 



<0,ad-i<Mpi. (54) 



It follows that 0^ ^ H'^Mp^, which confirms our approximation. Physically, this can also be 
understood from (46), describing a damped oscillator with potential energy pm0^/2M^. Once 
the field starts oscillating, the friction term can be ignored over the course of a few oscillations. 
The virial theorem then implies (0^) = pin(0^)/^^ < Pm0?ad-i/^^ ^ Pm, where 0rad-i ^ M by 
assumption. The symmetron kinetic energy is therefore subdominant to the matter component. 
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3.3 Matter-Dominated Era 

In the matter-dominated era, pm — SH'^Mp^, the scalar equation of motion becomes 

d^(J) 3d0 3M|, ^ 

+ oTi7 + ^r7F</' = 0, (55) 



dA^2 ' 2dA^ M2 

where A^ = ln(a/aeq) is the number of e-folds since matter-radiation equahty. Up to an irrelevant 
phase, the solution with initial condition = 0eq at equality is 



Again the field undergoes damped oscillations. By the onset of the phase transition {z ~ 1), the 
amplitude has dropped by 

Jt^-(r^-T-io-'. (57) 

|0eq| Vl + ^eq/ 

Interestingly, note from (52) and (56) that (p ~ a^^/^ both in the radiation- and matter- 
dominated eras. Combining (53) and (57), the total decay in the field amplitude from the moment 
it starts to evolve deep in the radiation-dominated era to the onset of the phase transition is 

J^.10-m"\ ,58) 

These results have been confirmed numerically. Figure 2a) shows the results of integrating 
the combined scalar field equation (40) and Friedmann equation (43) from z = 10^^ to z = 1, with 
0rad-i = 10~^M and M = 10~^Mpi. For this purpose, we approximated A{(f)) as the quadratic 
form (44) and chose fi = MpiHo/M and A = M^^H^/M^, as in (11) and (14), respectively, 
even though V{(j)) is irrelevant at these redshifts. As expected, the field remains stuck at its 
initial value until -Zunfreeze ~ 10^^, in agreement with (48). The field subsequently undergoes 
damped oscillations, with an amplitude decaying as (1 + 2;)^/^. By matter-radiation equality, the 
field amplitude has decreased by a factor of ~ 10^^, consistent with (53). During the matter- 
dominated era, the amplitude decays by an additional factor of ~ 10~^, which agrees with (57). 
Figure 2b) zooms in on the boxed region in Figure 2a), highlighting the fact that the oscillations 
are completely smooth, despite appearing jagged in Figure 2a). 
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Figure 2: Numerical solution of the symmmetron evolution as a function of redshift through 
the radiation- and matter-dominated eras, starting from 0rad-i = 10~^Mpi, with M = 10~'*Mpi. 
Figure (a) confirms that the field is initially frozen because of Hubble friction, but eventually 
undergoes damped oscillations. See main text for details. Figure (b) is a zoom-in on the boxed 
region, highlighting that the oscillations are smooth. 



4 Phase Transition and Late-Time Evolution 



In the previous Section, we have seen that starting from 0rad-i ^ M deep in the radiation- 
dominated phase, the symmetron subsequently evolves towards the symmetry-restoring point, 
culminating to a value ~ lO^'^(M/Mpi)^/^0i.ad-i by the onset of the phase transition. 

In this Section, we describe the symmetron dynamics during this late-time phase transition. 
For the quartic potential (1), as already mentioned in (15), the symmetron potential energy 
is insufficient to drive cosmic acceleration. Hence we will include a cosmological constant to 
reproduce the ACDM expansion history. Our main interest will be to check that the symmetron 
ends up near the symmetry-breaking value by the present time, as assumed in the analysis of 
spherically-symmetric solutions and tests of gravity. In Section 4.2 we will extend the analysis to 
potentials with greater promise of impacting both the recent expansion history and the growth 
history on linear scales today. 
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Figure 3: Evolution of the symmetron around the phase transition, a) At the onset of the 
transition (-Ztran = 1), the symmetron rapidly evolves towards (po, the minimum of the potential, 
and undergoes damped oscillations around this point. The dashed line indicates 4>mm, where we 
have chosen M = 10~^Mpi. b) Zoom-in on the boxed region, which shows that the oscillations are 
smooth. The solid line indicates (f)mm{z), the minimum of the effective potential given by (62). 



4.1 Quart ic Potential 



For this analysis, we add a constant to the quartic potential 1, 



(59) 



such that the minima have zero potential energy. As shown in (13), the minimum lies at 0o 
M^/Mpi <^ M, hence we can still approximate the couphng function as quadratic: ^(0) 
l + 0^/2M^. We also include a cosmological constant, with energy density p\ = SH^Mp^^l — Q 
where f2m denotes as usual the present fractional energy density in matter. The scalar equation 
of motion (40) then takes the form 



(0)n 
m ), 



+ 3ff0+^^^(4(l-nL°O+fiW(l + -f)0-/^^0 + A0^ 







(60) 



The phase transition occurs at a redshift Zt^an when the effective mass vanishes, that is. 



(61) 
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Figure 4: The fractional energy density in the symnietron for the quartic potentiaL 



For instance, the transition will happen at z^,. 
substituted Q 



1 for /x = y/l5HoMpi/M, where we have 



(0) 



0.25. This is consistent with fll). 



Since fi ^ Hq, shortly after the mass becomes tachyonic we expect the evolution of the 
symmetron to be adiabatic, in the sense that the field should track the minimum of the effective 
potential: 



0min(2;) 



/i 



3H^M^^ 



M2 



3^(0) W^i 
AM 



4 (1 - fiS^) + fiL°Hl + ^) 

1/2 



1/2 



:i + Ztr 



[l + z] 



(62) 



This "tracking" solution, valid for z < -Ztran, is an attractor because the effective potential is 
everywhere convex around this point. 

To check these results, we have numerically integrated the combined equation (40) and (43) 
from z = 3 until today, including a matter component and a cosmological constant, with fim = 
0.25. For the quartic potential, we chose /i = y/lbHoMpi/M, corresponding to Ztmn = 1, and 
A = Mp^Hq/M"^, consistent with (14). As initial conditions for the symmetron, we choose (j){z = 
3) = 10^^^ M, which matches the example of Figure 2. 

Figure 3 shows the evolution of the symmetron field as a function of redshift. At the onset 
of the transition, the symmetron rapidly rolls towards the minimum 0o = M^/Mpi (denoted by 
the dashed line) and oscillates around this point. Figure 3b) zooms in on the oscillations in the 
boxed region, confirming that the symmetron nicely tracks 0min(-2^); the minimum of the effective 
potential, as given by (62). 
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Figure 5: The generalized potential studied in Section 4.2, shown here for illustrative purposes 
with M = O.lMpi. The minima occur at ~ ±M. By construction, the height of the potential is 
~ i^gMpj, and the mass around the minima is ~ Hq. 

Figure 4 shows the fractional energy density in the symmetron. This confirms that the 
symmetron backreaction is negligible, as assumed in the analytical approximations above. The 
resulting expansion history is therefore indistinguishable from ACDM in the quartic case. 



4.2 Generalized Potential and Coupling Function 

To maximize the impact on the late-time expansion and growth histories, we generalize V{(f)) such 

that the height of the potential is comparable to the dark energy scale, Vq ~ H^Mp^, and the mass 

around the minimum is of order Hubble, mo ~ Hq. The following general form satisfies these 

requirements: 

/ i\/r \ 

(63) 



y(0) = H^M^, e 



2A/f2 -acf>y2M^ 



+ 



M 



>2/2M2 



M, 



PI 



where a is a constant. This is sketched in Figure 5. For |0| <^ M, 



Vi<t>)^H^oMi,-l^^^<^' + 



(64) 



2 M2 

this matches (1) to quadratic order, up to a constant, with /i = y/aHoMpi/M, which is consistent 
with (11). 

The minima in this case are located at values of order M , 

„i/ 



'6 ,,, 1/2 /«^/^Mpi 
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(65) 



Although one would have to carefully revisit the various tests of gravity analyzed in Section 2.3, 
a naive application of the thin shell condition ARq/Rq = 0o/'^g^pi ^ 1 for our galaxy would 
imply a tighter bound M < 10~^Mpi. We will keep this bound in mind, though our conclusions 
will be insensitive to the precise choice of M. 

Note that the potential evaluated at the minima is V{(f)o) — H^MMpi ^ H^Mp^i hence the 
difference in potential energy with the maximum, AV ~ H^Mp^, is of the order of the dark energy 
scale, as desired. Meanwhile, the mass of small fluctuations around the minima is 



With M = 10~^Mpi and a = 1, this gives mo ~ 10~^Ho, which is smaller than Hubble. Hence, as 
desired, the symmetron-mediated force extends to the largest observable scales. Thus, as Figure 5 
illustrates, the potential has a sharply-peaked maximum at = and shallow minima at = ±0o- 

Because reaches values of order M in this case, the quadratic approximation (44) for A{(j)) 
breaks down. Not only are corrections expected for (p ~ M, they are in fact necessary to prevent 
the symmetron-mediated force from being unacceptably strong. Indeed, if A{(j)) = 1 + 0^/2M^ 
were valid full-stop, from (12) the force between two test masses in vacuum would be F^/F-^ ~ 
Mpi/M"^ ~ 10^^. Hence we must generalize A{(f)) by demanding that the scalar force be at most 
comparable to gravity for all field values: d\nA/d(f) < 1/Mpi for all 0. As a particular example, 
consider 

This reduces to ^4(0) ~ 1 + 0^/2M^ for small (f) and tends to ^4(0) ~ 1 + |0|/Mpi for large 0. 

With this choice of V{(f)) and A{(f)), we numerically solved the cosmo logical equations (40) 
and (43). Since our goal is to have the symmetron drive cosmic acceleration, we did not included a 
cosmological constant in this case. To facilitate the comparison with the quartic case of Section 4.1, 
we chose M = 10~'^Mpi and picked a so that the transition occurs at Ztra,n — 1- 

The solution for the scalar field is shown in Figure 6. Shortly after the onset of the phase 
transition, we notice that the symmetron overshoots the minimum and reaches a value of order 
Mpi. Consequently, the field has not yet converged to 0o by the present time — note that the 
Figure even extends to ^ < — which is problematic for tests of gravity. 

Figure 7 shows the Hubble parameter H{z) for the solution (solid line). For comparison, 
the dotted line is a ACDM model with ilm = 0.25, while the dashed line is an Einstein-de Sitter 
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Figure 6: The behavior of after the transition for the generahzed potential (63) and couphng 
function (67). The dashed hne denotes the minimum (pQ, given by (65), where we have taken 
M = 10~^Mpi. We have chosen parameters such that ^tran = 1- Note that the plot extends 
to future redshifts [z < 0). Because of the asymmetric form of the potential, in this case the 
symmetron has not yet settled to the minimum by the present time. 



(r^m = 1) universe. Aside from the small blip at 2; ~ ^tran, the solution rapidly converges to 
the Einstein-de Sitter evolution, indicating that the symmetron energy density is negligible for 
z < Ztran- The failure to drive cosmic acceleration is due to the sharply peaked maximum, which 
leads to a rapid completion of the transition and the field overshooting the minimum. 

An obvious way to avoid the overshooting problem requires tuning parameters such that the 
transition occurs very recently (2;tran <^ !)• In this way, one can make the field value today of order 
M, as the field evolves towards much larger values. At the same time, the potential energy at 
the maximum can drive cosmic acceleration until the transition, matching the ACDM expansion 
history until the very recent past. This solution, though consistent, is clearly not appealing as it 
requires delicate tuning of parameters. To summarize, though we have focused on a particular class 
of potentials, it appears difficult for the symmetron to drive cosmic acceleration while satisfying 
local tests of gravity, unless one considers tuning parameters. 
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Figure 7: Hubble parameter as a function of redshift for the symmetron with generahzed po- 
tential (63) and coupling function (67) (solid line), where Hq = 72 knis~^Mpc~^ denotes the 
measured value. For comparison, the thick dashed line is a ACDM cosmology with fim = 0.25, 
while the thin dashed line is an Einstein de Sitter {ilm = 1) universe. In all three cases, we have 
fixed the expansion history deep in the matter-dominated era. 



5 The Inflationary Phase 



Having studied the symmetron evolution throughout the radiation- and matter-dominated era, 
including the late-time symmetry-breaking transition, we now go back and study more closely the 
issue of initial conditions deep in the radiation era. While the initial value 0rad-i at the onset 
of the radiation-domination phase was treated as a free parameter in Section 3, in this Section 
we will see how this initial condition comes about from the field dynamics during inflation. We 
will find that the field amplitude decays by an exponential amount during inflation, but can be 
exponentially resurrected during reheating. In particular, the condition 0rad-i ^ M assumed in 
Section 3 will translate to a constraint on the reheating phase. 

Of course, the inflationary analysis below must itself assume an initial value 0inf-i, which 
depends on the pre-inflationary history. To be conservative, however, we will set 0inf-i — M, 
corresponding to the largest field value consistent with the symmetron energy density being sub- 
dominant during inflation. 
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5.1 Dynamics During Inflation 

For simplicity, we treat the inflationary phase as exact de Sitter expansion, a{t) ~ e^'^'"'*, neglecting 
the backreaction of the symmetron field. By virtue of its universal coupling to all matter fields, 
the symmetron couples to the infiaton, being sourced by Tinf = (—1 + 3i(;inf)pinf — — 4pinf- 

The assumption that the symmetron is a spectator is of course not necessary — though 
its backreaction could in principle prevent inflation from starting — but is made to simplify the 
analysis. In this case, as in the discussion of the radiation-dominated evolution, for consistency 
the initial field value 0inf-i at the onset of inflation must satisfy 0inf-i ^ M, such that (44) applies. 
In particular, to a good approximation A{(f)) ~ 1, and the infiaton contribution to the Friedmann 
equation (43) is nearly 0-independent and hence constant: 

3Hl,Ml, ^ pinf . (68) 

Meanwhile, the scalar equation of motion (40) reduces to 

^ + 3i7inf0 + 12^^^^^0 = . (69) 

Again, this is just the equation for a damped harmonic oscillator, with solution 

= 0inf-ie-=^^-f*/'cos(u;t) , (70) 

where 

, = 2y^^!s^Xn7Zy^2V3^^ei^. (71) 

M V 48 \MpJ M ^ ^ 

Clearly, u ^ -ffinf , and the symmetron oscillates rapidly as inflation proceeds. The choice of phase 
in (70) corresponds to the symmetron starting nearly from rest from 0inf-i at the onset of inflation 
(t = 0). 

Since the effective symmetron mass is large compared to Hubble, its energy density averaged 
over many oscillations redshifts as dust: 

(p.) - l^'<l>L-.e-'"-^' = &Hl,Ml, (^) ' e-^^-* . (72) 

This confirms that 0inf-i ^ M is indeed necessary to neglect the symmetron backreaction from 
the outset, an approximation which only gets better in time, as (72) implies. Over the entire 
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period of inflation, therefore, the fleld amphtude decays by 

I'/'inf-il 

where N^^f is the number of inflationary e- folds. For the standard value A^i^f = 60, corresponding 
to GUT-scale inflation, this gives a suppression factor of e"^*^ ~ 10~^^. 

5.2 Reheating 

We next track the evolution of the symmetron during reheating, as the inflaton decays into radi- 
ation. Once again neglecting the symmetron backreaction, the scalar equation (40) becomes 

4> + 3H<P+^{l-3w)<l) = 0, (74) 

where we have also neglected the potential term V^. Here p and w respectively denote the energy 
density and equation of state of the combined inflaton-radiation fluid. To simplify the analysis, 
we will assume that reheating occurs rapidly compared to a Hubble time, so that H ~ Hint 
remains approximately constant. Moreover, in this regime the Hubble friction term can be ignored, 
and (74) reduces to 

0+^(1 -3w7)0 = O, (75) 

where we have used p = SH'^Mp^ ~ SH^^^Mp^, and where u ~ 2^/3MplHinl/M was deflned in (71). 

To proceed, we need to specify w{t) during reheating. For concreteness, we consider the 
simplest case of the inflaton decaying perturbatively into radiation as it oscillates around the 
minimum of its potential. In the process, w oscillates between —1 and +1, and approaches 1/3 as 
the universe becomes dominated by radiation. To model this evolution, consider 

w{t) = ^ [1 - e-^'"f* (1 + 3 cosminft)] , (76) 

where t = now marks the onset of the reheating phase. This form is such that the equation of 
state starts out with w{0) = —1, thereby matching smoothly to the inflationary phase. At late 
times, w — )■ 1/3, as desired. Manifestly, Fjnf and minf denote the inflaton decay rate and mass, 
respectively. 

With this choice of w, (75) becomes 

4> + — e-^-f* (1 + 3 cos minft) 4> = . (77) 
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Figure 8: Symmetron evolution during reheating for different parameter values. In each case, 
the thick line indicates the numerical solution, whereas the thin line is the Mathieu cosine func- 
tion (80). All plots have uj = lO'^Hi^f and M = 10~'*Mpi. a) In this case, Finf ^ fn-u^i, hence 
reheating has essentially completed before the tachyonic amplification can kick in. b) In the op- 
posite regime, Fjnf ^ mi^f, reheating is slow. Since minf = u lies in the unstable region of the 
Mathieu cosine solution, the field is exponentially amplified. The amplification is eventually cur- 
tailed once the inflaton decays and reheating takes place. The dashed line indicates the maximum 
allowed amplification, corresponding to 0rad-i — M. c) This case also has F^f ^ m-mf, as in b), 
but TTiinf = lO^ci; lies in the stability region of the Mathieu cosine. Note the logarithmic scale on 
the time axis in this case. 
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Thus the evolution is characterized by three dimensionful parameters: the inflaton decay rate r^f , 
the inflaton mass rrii^f, and the effective symmetron mass u during inflation. Since reheating is 
assumed to proceed rapidly on a Hubble time, all three parameters are ^ -f^inf- Whenever the 
equation of state achieves w > 1/3, the effective symmetron mass term in (77) becomes tachyonic, 
which can amplify the field value. As w -^ 1/3, the effective mass is driven to zero, and the field 
evolves freely. Eventually the field will be slowed down by Hubble friction. 

The behavior of the solution can be understood by first considering the regime t <C Frj:, 
where the exponential damping term is nearly constant. Defining x = minft/2 and y = (p/(pmi-end, 
this scalar equation reduces to a Mathieu equation 

d'y 



dx"^ 
with 



+ [a-2qcos2x]y = 0, {7i 



LO^ 



a=^] q = --a. (79) 

The unique solution with y{0) = 1 and dy/dx{0) = is the Mathieu cosine function, y = C{a, q, x), 
which in our case is given by 

(Pit) ~ 0i,f_end C (a, -^a, ^"j for t « TrJ . (80) 

The stability of the symmetron evolution is therefore determined by the stability properties of 
the Mathieu cosine function. For a < 10"^, corresponding to mjnf > 10 w, we find that the solution 
undergoes stable oscillations, consistent with the fact that C{a,q,x) ~ cos{y/ax) for |g| <^ 1. For 
a>10~^, corresponding to minf <10a;, the solution is unstable, and the symmetron amplitude 
grows exponentially. The growth is eventually cut off when t ~ Fj^^^, and the source becomes 
exponentially suppressed, indicating that reheating has nearly completed. 

Figure 8 confirms these expectations. In all cases, the thick line is the numerical solution, 
while the thin line is the Mathieu cosine function (80). Figure 8a) is an example of rapid reheating, 
Finf ^ TTiinf- In this case, reheating is essentially complete before the tachyonic amplification can 
kick in, and the Mathieu cosine is a good approximation for less than one symmetron oscillation. 

Figures 8b) and c) are examples of slow reheating, F^f ^ ^Tijnf, in which case the Mathieu 
cosine matches the exact solution for many oscillations. Figure 8b) corresponds to r/iinf = u, which 
lies within the unstable region of the Mathieu cosine parameter space. The field is amplified by 
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many orders of magnitude, until its growth is tamed after a time t ~ Fj^J when reheating takes 
place. The dashed line indicates the maximum allowed amplification so that |0rad-i| < M, as 
discussed shortly. The parameters have been chosen in this case such that the amplification 
barely satisfies this bound. Figure 8c) corresponds to minf = lO^ci;, which lies within the stable 
region of the Mathieu cosine. Note that the scale on the t-axis is logarithmic in this case. 

After several decay times, Finft ^ 1, the driving term becomes negligible, and (77) reduces 
to </) ~ 0. Thus the field eventually grows linearly with time, as confirmed in all examples of 
Figure 8. Of course this linear growth eventually gets cut off by Hubble friction. The asymptotic 
value is is identified with 0rad-i, the initial value for the radiation-dominated evolution. 

We can derive a constraint on the reheating parameters by demanding that 0rad-i < M, 
for consistency of neglecting the symmetron backreaction, as assumed throughout Section 3.2. 
For simplicity, we estimate 0rad-i as the field value at t ~ Fj^j, and ignore the relatively small 
subsequent growth. To be most conservative, we take 0inf-i = M at the onset of inflation. 
Assuming A^jnf = 60 e-folds of inflation, (73) imphes 0inf-end — 10~^^M. The condition 0rad-i < M 
therefore corresponds to 



<f>{t = rrJ) 



< 10^^ . (811 

'inf— end 

As discussed above, there is minimal growth if Finf ^ minf (in which case reheating happens 
quickly) and/or mjnf >10co (in which case the Mathieu cosine is stable). In the regime Fjnf ^ rnjnf 
and TTiinf < 10 w, however, we have found numerically that (81) is satisfied for 

N 0.7 

^) . (82) 

mini J 

To summarize, due to possible parametric resonance during reheating, the symmetron amplitude 
can increase exponentially. Demanding that this amplification remains within bounds imposes 
constraints on the parameters, which we have derived in the limit that all time scales involved are 
short relative to a Hubble time. 

6 Conclusions 

Despite the overwhelming evidence for the existence of dark energy and dark matter, little is known 
about their underlying fundamental physics. Over the last few years there has been considerable 
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activity exploring the possibility that the dark sector includes new light degrees of freedom (gen- 
erally scalar fields) that couple not only to dark matter but also to ordinary (baryonic) matter. 
Much of the research efforts have focused on the development of screening mechanisms to explain 
why such scalars, if light, have escaped detection from laboratory/solar system tests of gravity. 
The manifestation of these scalar fields therefore depends sensitively on their environment, which 
in turn leads to striking experimental signatures. 

One such screening mechanism is the symmetron mechanism, which relies on a scalar field 
having a VEV that depends on the ambient matter density. The VEV is small in regions of high 
density, and large in regions of low density. In addition, the coupling of the scalar to matter is 
proportional to this VEV, hence the scalar couples with gravitational strength in regions of low 
density, but couples much more weakly in regions of high density. 

In this paper we have derived the cosmological expansion history in the presence of a sym- 
metron field, tracking the evolution of the scalar through the inflationary epoch, the phase of re- 
heating, the standard radiation- and matter-dominated eras, and through the late-time symmetry- 
breaking phase transition. For a wide range of initial field values at the onset of inflation, we have 
shown that the symmetron ends up in the symmetry-breaking vacuum by the present time, as 
assumed in the derivation of static, spherically-symmetric solutions and analysis of tests of gravity. 

For the fiducial quartic potential, the scales involved are too small to drive late-time cosmic 
acceleration, hence we included a cosmological constant. The symmetron backreaction is consis- 
tently small throughout the evolution, hence the expansion history is consistent with observations. 
We also introduced a general class of potentials, engineered so that the potential energy difference 
is of order the present critical density, and the mass around the minimum is of order Hubble today. 
In the simplest examples considered here, however, the symmetron was unsuccessful in driving 
cosmic acceleration. 

Our results form the foundation of future investigations of the cosmological implications 
of symmetrons. The natural next step is to study the evolution of density perturbations in 
the presence of a symmetron, both in terms of the linear growth history and in the non-linear 
regime using N-body simulations. As in the chameleon and Vainshtein screening mechanisms, 
the additional scalar force should enhance structure growth at late times. It would be interesting 
to contrast the impact on structure formation among different screening mechanisms and seek 
distinguishing signatures. 
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Another cosmological probe of symmetron physics is the production of topological domain 
walls. In this paper we have circumvented this fascinating issue by implicitly focusing on initial 
conditions where the scalar field ends up in the vicinity of = at late times, but sufficiently 
displaced such that the phase transition can proceed homogeneously. More generally, we expect 
the symmetron to result in late-forming domain wall networks. It will be interesting to study their 
formation and contrast the phenomenology of symmetron defects with standard topological defect 
networks by comparing their observational impact, e.g., on CMB anisotropy, structure growth and 
gravitational lensing. 
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